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Abstract 

We study analytically the initial value problem for a self-interacting (mas- 
sive) scalar-field on a Reissner-Nordstrom spacetime. Following the no-hair 
theorem we examine the dynamical physical mechanism by which the self- 
interacting (SI) hair decays. We show that the intermediate asymptotic be- 
haviour of SI perturbations is dominated by an oscillatory inverse power-law 
decaying tail. We show that at late-times the decay of a SI hair is slower than 
any power-law. We confirm our analytical results by numerical simulations. 

I. INTRODUCTION 

The late-time evolution of various fields outside a collapsing star has an important im- 
plications to two major aspects of black-hole physics: the no-hair theorem and the mass- 
Inflation scenario. The no-hair theorem, introduced by Wheeler in the early 1970s, states 
that the external field of a black-hole relaxes to a Kerr-Newman field characterized solely by 
the black-hole's mass, charge and angular-momentum. Thus, it is of interest to explore the 
dynamical physical mechanism responsible for the relaxation of perturbations fields outside 
a black-hole and to determine the decay- rates of the various perturbations (which differ from 
one field to the other). The mechanism by which massless neutral fields are radiated away 
was first studied by Price fl[]. The physical mechanism by which a massless charged scalar 



hair is radiated away was studied in @,[|. In this paper we study the physical mechanism 
responsible for the decay of a self-interacting (massive) scalar-hair. 

The asymptotic late-time tails along the outer horizon of a rotating or a charged black- 
hole are used as initial input for the ingoing perturbations which penetrates into the black- 
hole. These perturbations are the physical cause for the well-known phenomena of mass- 
inflation [§J. In this context, one should take into account the existence of massive tails out- 
side the collapsing star. Here we study analytically the intermediate asymptotic behaviour 
of such self-interacting (massive) perturbations fields. We study the late-time asymptotic 
behaviour numerically and confirm the numerical results of Burko || . 

The plan of the paper is as follows. In Sec. FT] we describe the physical system and 
formulate the evolution equation considered. In Sec. |TTT] we formulate the problem in terms of 
the black-hole Green's function using the technique of spectral decomposition (this section is 
analogous to the one given in 0). In Sec. |TV| we study the intermediate asymptotic evolution 
of SI scalar perturbations on a Reissner-Nordstrom background. We find an oscillatory 
inverse power-law behaviour of the perturbations at a fixed radius and along the black-hole 
outer-horizon. We find that the dumping exponents which describe the intermediate fall-off 
of SI perturbations are smaller compared with the massless (neutral) dumping exponents. 
In Sec. |V| we verify our analytical results by numerical simulations. We conclude in Sec. |V| 
with a brief summary of our results and their implications. 



II. DESCRIPTION OF THE SYSTEM 

We consider the evolution of SI scalar perturbation fields outside a collapsing star. The 
external gravitational field of a spherically symmetric collapsing star of mass M and charge 
Q is given by the Reissner-Nordstrom metric 



, 2M Q' 2 \ , , / 2M Q 2 \ 1 , „ 
ds 2 = - 1 + ^ T )dt 2 +[l + dr 2 + r 2 dVl 2 . (1) 
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Using the tortoise radial coordinate y, defined by dy = dr/X 2 where A 2 = 1 — — + ^r, the 
metric becomes 



ds z = X\-df + dy')+r 2 dn 2 



(2) 



The wave equation for the SI scalar-field is 



w - 1/ 



o , 



(3) 



where ?7(|0| 2 ) is the self-interaction potential and £7'(|0| 2 ) = <i[/(|0| 2 )/<i0*. Since we study 
the evolution of small perturbations we approximate U'(\<f)\ 2 ) by m 2 (we assume that m 
is real), neglecting terms of higher order in cf). Resolving the field into spherical harmonics 



= J2 r ) Y™^, V 9 )/ 7 " one obtains a wave-equation for each multiple moment 

l,m 
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where 



V = V M ,Q,l,m (r) = 1 
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III. FORMALISM 



The time-evolution of a SI scalar-field described by Eq. is given by 

^(z/> *) = / x ; *)^tfa 0) + G t(y, x; t)ip(x, 0)] dx , 
for t > 0, where the (retarded) Green's function G(y, x\ t) is defined as 

d 2 d 2 



V(r) 



G(y,x;t) = 5(t)5(y-x) 



(6) 



(7) 



dt 2 dy 2 

The causality condition gives us the initial condition G(y,x;t) = for t < 0. In order to 
find G(y, x; t) we use the Fourier transform 

n 

(8) 



G(y,x;w) = / G(y, x;t)e lwt dt . 
Jo- 

The Fourier transform is analytic in the upper half w-plane and it satisfies 

[if* +n}2 G (v> x;w) = 6(y-x) . 



(9) 



G(y, x; t) itself is given by the inversion formula 

"| roc+ic _ 



G(y, x;t) = — G(y, x; w)e~ lwt dw , (10) 

ZTT J — oo+ic 

where c is some positive constant. 

Next, we define two auxiliary functions ip\(y,w) and ijj2{y,w) which are linearly inde- 
pendent solutions to the homogeneous equation 

(J^+u?-v\j>i(y,w) = 0, i = l,2. (11) 

Let the Wronskian be 

W(W) = W$1, fc) = Tpl<4>2,y - ^l,y , (12) 

where W(w) is y- independent. Using the two solutions ipi and ip 2 , the black- hole Green's 
function can be expressed as 

N 1 \ $i(y,w)<ip2(x,w) , y<x, 

G(y,x-w) = -——( _ _ (13) 

1 > y iJi(x,w)ip 2 (y,w) , y > x . 

In order to calculate G(y,x;t) using Eq. fllPf), one may close the contour of integra- 
tion into the lower half of the complex frequency plane. Then, one finds three distinct 
contributions to G(y,x;t) [|] : 

1. Prompt contribution. This arises from the integral along the large semi-circle. It is this 
part, denoted G F , which propagates the /wg/i-frequency response. For large frequencies 
the Green's function approaches to the one of a massless scalar-field on a flat spacetime 
background. This term contributes to the short-time response and can be shown to 
be effectively zero beyond a certain time. Thus, it is not relevant for the late-time 
behaviour of the field. 

2. Quasinormal modes. These arise from the distinct singularities of G(y,x;w) in the 
lower half of the complex w-plane and is denoted by G®. These singularities occur at 
frequencies for which the Wronskian (|l^) vanishes. G® is just the sum of the residues 



at the poles of G(y, x; w). Since each mode has Imw < it decays exponentially with 
time. 

3. Tail contribution. As will be shown later the intermediate asymptotic tail is associated 
with the existence of a branch cut (in ip 2 ) placed along the interval —m < w < m. 
This tail arises from the integral of G(y, x; w) around the branch cut (denoted by G ) 
which leads to an oscillatory inverse power-law behaviour of the field. Since we are 
interested in the intermediate asymptotic behaviour of a SI scalar-field our goal is to 
evaluate G°\y,x;t). 

IV. THE INTERMEDIATE ASYMPTOTIC BEHAVIOUR OF A 
SELF-INTERACTING SCALAR-FIELD 

A. The M < r < M/(Mm) 2 approximation 

It is well known that the late-time behaviour of massless fields is determined by the 
backscattering from asymptotically far regions 0,0,0 . Thus, the late-time tails of massless 
fields are dominated by the /otofrequencies contribution to the Green's function, for only 
low frequencies will be backscattered by the small spacetime curvature or by the small 
electromagnetic interaction at these asymptotic regions. On the other hand, it is also well 
known that massive-tails exist even in a flat spacetime ||. This phenomena is related to the 
fact that different frequencies forming a massive wave packet have different phase velocities. 
As will be shown in this paper, the intermediate asymptotic behaviour of a SI scalar-field 
on a Reissner-Nordstrom background is dominated by flat spacetime effects. Namely, at 
intermediate times the backscattering from asymptotically far regions (which dominates the 
tails of massless fields) is negligible compared to the flat spacetime massive tails that appear 
here. 

Let us assume that both the observer and the initial data are situated far away from 
the black-hole. We expend the wave-equation ([□]) for the SI scalar-field (in the field of 
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the black- hole) as a power series in M/r and Q/r and obtain (neglecting terms of order 
0[(Mm)2] an d higher) 

d 2 



dr' 



+ w' 



9 AMw 2 - 2Mm 2 1(1 + 1) 
m + 



(14) 



where £ = A^>. The term proportional to M/r represent the Newtonian potential. If we 
further assume that both the observer and the initial data are situated in the region r <C 
M/(Mm) 2 (and M<r), and we are interested in the intermediate asymptotic behaviour 



of the field (r <C t <C M / (Mm) ) we can further approximate Eq. (|14j ) by 



dr 1 



/(/ + !) 



(15) 



Replacing Eq. (|T4|) with Eq. dl5|) means that we neglect the backscattering of the field 
from asymptotically far regions. Thus, the intermediate asymptotic behaviour of SI scalar 
perturbations on a black-hole (or a star) background depends only on the field's parameters 
(namely, on the mass of the field) and it does not depend on the spacetime parameters. The 
validity of this conclusion is verified by numerical simulations (see Sec. |V|). 
Let us now introduce a second auxiliary field <fi defined by 



£ = r l+1 e Ww2 ~ m2r 4>(z) 



(16) 



where 



z = —2iwr . 



(17) 



(z) satisfies the confluent hypergeometric equation 



d 2 



d 



z— + (21 + 2- z)— -{1 + 1) 

dz z dz 



:z) = o. 



(18) 



The two basic solutions required in order to build the black-hole Green's function are (for 
\w\ < m) 



$ x = Ar l+1 e~ mr M(l + 1,2/ + 2, 2wr) 



(19) 



and 



^ 2 = Br l+1 e-™ r U(l + 1,2/ + 2, 2wr) 



(20) 



where w = \Jm 2 — w 2 . A and B are normalization constants. M(a, b, z) and U(a, b, z) 
are the two standard solutions to the confluent hypergeometric equation ||. U(a,b,z) is a 
many-valued function, i.e. there is a cut in Using Eqs. 13.6.6 and 13.6.24 of || one may 
write these solutions in a more familiar form 



and 



A = ^r(/ + 2)(i 07 )-(H-i)j I+J ( ro , 



^ 2 = n-^Br^(2zu)- {l+ ^K l+ i(wr) 



(21) 



(22) 



where and K l+ i are the modified Bessel functions. 



Using Eq. (10), one finds that the branch cut contribution to the Green's function is 
given by 



III J-m 



' iwt dw . 



(23) 



W{we™) W(w) 

For simplicity we assume that the initial data has a considerable support only for r-values 
which are smaller than the observer's location. This, of course, does not change the late-time, 
behaviour. 

Using Eqs. 9.6.30 and 9.6.31 of |§, one finds 



and 



0i(r, we m ) = ?/>i(r, w) 



B7r3(-l) , + 1 2- a ~ 
■02 (r, roe 7 ™) = -^ 2 (r, ro) + — — rn ; 3N ^i(r, cs) 



(24) 



a r(z + §) 

Using Eqs. ( p^) and (p3[) it is easy to see that 



(25) 



(26) 



From which we obtain the relation 



W{we™) W(w) A r(/ + §) W(w) ' 1 j 

Since W(w) is r- independent, we may use the wr — > asymptotic expansions of the Bessel 
functions (given by Eqs. 9.6.7 and 9.6.9 in ||) in order to evaluate it. One finds 

W{w) = -^ABn-^(2l + l)T(l + i)ro- (2m) . (28) 



Finally, substituting ( p7[) and ( |28|) in (^3|) we obtain 

(-1) Z+1 2 2 - 2Z 
A 2 (2/ + l)r(Z + |)r(/ + |) 



G C (y,x;t) = — - ,\ J„ ■ r — - ■ 3 , / My,w)Mx,w)w 2l+1 e- M dw . (29) 



B. Intermediate asymptotic at a fixed radius 

It is easy to verify that in the large t limit the effective contribution to the integral in 
arise from \w\ —0(m — |) or equivalently w = O(y^). This is due to the rapidly oscillating 
term e~ lwt which leads to a mutual cancelation between the positive and the negative parts 
of the integrand. In order to obtain the intermediate asymptotic behaviour of the field at a 
fixed radius (where x, y <C t), we may use the wr <C 1 limit of ipi(r,w). Using Eq. 9.6.7 
from H one finds 

V>i(r,si7) ~ ^Ar l+1 . (30) 

We obtain 

G c (y,x-t) = ^=7^^ ' (31) 

which in the t 3> m _1 limit becomes 

G C (y,x;t) = ^^m l+ Hxy) l+1 t-^cos[mt - {\l + ^tt] . (32) 

Thus, the intermediate asymptotic behaviour of the SI field at a fixed radius is dominated 
by an oscillatory inverse power-law tail. 
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C. Intermediate asymptotic along the black- hole outer horizon 



Next, we consider the behaviour of the SI scalar-field at the black- hole outer-horizon 



r + . While Eqs. ([H]) and (22) are approximate solutions to the wave-equation ([11]) in the 
M ■€ r < M/(Mm) 2 region, they do not represent the solution near the horizon. As 
y — > — oo the wave-equation ( O ) can be approximated by the equation 

i> m + w 2 i) = . (33) 

Thus, we chose 

^{y,w) = C{w)e- iw \ (34) 

and we use the form (|30|) for ip\(x, w). In order to match the y <C — M solution to the 
y ^> M solution we assume that the two solutions have the same temporal dependence. 
This assumption has been proven to be very successful for massless neutral Jl(| and charged 



Pll perturbations. In other words we assume that C(w) is w- independent. In this case one 
should replace the roles of x and y in Eqs. (|23|) and (p9|). Using fl29|), we obtain 



G c (y, x- 1) = T^y^m l ^y^v-^cos[mv - (h + ^tt] , (35) 

where T is a constant. 

Thus, the intermediate asymptotic behaviour of the SI field along the black-hole outer- 
horizon is dominated by an oscillatory inverse power-law tail. 



V. NUMERICAL RESULTS 

It is straightforward to integrate Eq. (f|) using the method described in []10[ . The late- 
time evolution of a SI scalar-field is independent of the form of the initial data used. The 
results presented here are for a Gaussian pulse on u = 

ip(u = 0, v) = Aexp j— [(v — v )/a} 2 } , (36) 



where the amplitude A is physically irrelevant due to the linearity of Eq. (f|). It should be 
noted that the evolution equation (^) is invariant under the rescaling 

r — > ar , t — > at , M — > aM , Q — > aQ , m — > m/a , (37) 

where a is some positive constant. The black- hole mass and charge are set equal to M = 0.5 
and Q = 0.45, respectively. We have chosen an initial field-profile with m = 0.01 ,vq = 50 
and <j = 2. We studied the behaviour of the field ip at a fixed radius and along the black-hole 
outer horizon (approximated by the null surface u = u max , where u max is the largest value of 
u on the grid). The numerical results for the I = mode are shown in Fig. [I]. Initially, the 
evolution is dominated by the prompt contribution and the quasinormal ringing. However, 
at intermediate times a definite oscillatory power-law fall off is manifest. The power-law 
exponents are —1.48 at a fixed radius (top curve) and —1.47 along the black-hole outer- 
horizon (bottom curve). These values are to be compared with the analytically predicted 
value of —1.5, see Eqs. (^) and (pop. The period of the oscillations (for |^|) equals 
T = ir/m to within 0.1%, again in agreement with the predicted value. Fig. 0. depicts the 
dependence of the intermediate asymptotic tails (at a fixed radius y = 50) on the multipole 
index /. The numerical values of the power-law exponents, describing the fall-off of the field 
at intermediate times (after a period of quasinormal ringing) are —1.49, —2.50, —3.50 and 
—4.51 for Z=0,l,2 and 3, respectively. These numerical values are in excellent agreement 
with the analytically predicted values of —1.5, —2.5, —3.5 and —4.5. Again, the period of 
the oscillations is T = ir/m to within 0.1%, in agreement with the predicted value. 

The analytical derivations and their numerical confirmations presented so far are re- 
stricted to the intermediate asymptotic regime. The late-time evolution of the field (for the 
/ = mode) is shown in Fig. |3] (for the clarity of the figure we display only the maximas of 
the oscillations). The initial data are those of Fig. |1|. Shown are the behaviour of the field 
l^l along the asymptotic regions of timelike infinity i + (top curve) and along the black-hole 
outer horizon (bottom curve). It is clear that the field's amplitude decays with time, in 
agreement with the no-hair theorem. The decay rate is slower than any power-law. Again, 
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we find that the field oscillates with a period of T = ir/m to within 0.2%. We have found 
that the larger is the field's mass the sooner it leaves the intermediate asymptotic phase of 
an inverse power-law decay. 



VI. SUMMARY AND PHYSICAL IMPLICATIONS 

We have studied analytically the intermediate asymptotic evolution of a SI scalar-field on 
a Reissner- Nordstrom background. Following the no-hair theorem we have focused attention 
on the physical mechanism by which a SI hair decays. The main results and their physical 
implications are: 

Oscillatory inverse power-law tails develop at intermediate times at a fixed radius and 
along the black-hole outer horizon (as long as the initial data has a considerable support 
only in the region M«rC (Mm) 2 - Actually, our analytical derivations hold even in the 
case M <^ r <^ ^Mmy> ^ or ^ > 0)- The dumping exponents, describing the fall-off of a ST field 
at intermediate times, are smaller compared with those of massless neutral perturbations 
fiyiOfl. For / > 77j|eQ| — \ (where e is the field's charge) these dumping exponents are also 



smaller than those of massless charged perturbations ||. While the asymptotic behaviour of 
massless perturbations is dominated by backscattering from asymptotically far regions [and 
thus depends on the spacetime parameters, M (for neutral perturbations) and Q (for charged 
perturbations)], the intermediate asymptotic behaviour of a SI field depends on the field's 
parameters (namely, on the field's mass m). In other words, dealing with SI perturbations, 
one may neglect the backscattering from asymptotically far regions at intermediate times. 

Using the numerical scheme we have shown that at late-times the inverse power-law 
decay is replaced by another pattern of decay, which is slower than any power-law. This 
late-time behaviour deserves a further analytic study. The late-time behaviour of SI per- 
turbations implies that a black-hole which forms from a gravitational collapse of SI fields 
becomes "bald" slower than one which forms during a gravitational collapse of massless 
fields. Since SI perturbation fields decay at late-times slower than any power-law they are 
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expected to cause a mass-inflation singularity during a gravitational collapse which leads 
to the formation of a black-hole. Moreover, SI fields are expected to dominate the mass- 
inflation phenomena during a gravitational collapse (this is caused by the slower decay of 
SI perturbations compared with massless ones). 

We have confirmed our analytical results by numerical calculations, showing that the 
intermediate asymptotic behaviour of a SI field is dominated by an oscillatory inverse power- 
law tail (with the analytically predicted dumping exponents and oscillation frequency). 
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FIG. 1. Evolution of the SI field on a Reissner-Nordstrom background, for 

I = 0, M = 0.5, Q = 0.45 and m = 0.01. The initial data is a Gaussian distribution with vo = 50 
and a = 2. The field at a fixed radius (y = 50) is shown as a function of t. Along the black-hole 
outer horizon the field is shown as a function of v. The oscillatory power-law fall off is manifest at 
intermediate times. The power-law exponents are —1.48 at a fixed radius (top curve) and —1.47 
along the black-hole outer- horizon (bottom curve). These values are to be compared with the 
analytically predicted value of —1.5. The period of the oscillations is T = ir/m to within 0.1%, in 
agreement with the predicted value. 
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FIG. 2. The amplitude of the field \ip{y = 50,i)| for different multipoles Z=0,l,2 and 3 
(from top to bottom). The power-law exponents are —1.49,-2.50,-3.50 and —4.51, respectively, 
in excellent agreement with the analytically predicted values of —1.5,-2.5,-3.5 and —4.5. The 
oscillations period is T = ir/m to within 0.1%, in agreement with the predicted value. The initial 
data are those of Fig. ||. 
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FIG. 3. Late-time evolution of the SI field \ip\ on a Reissner-Nordstrom background. The field 
at future timelike infinity (y = 50) is shown as a function of t (top curve). Along the black-hole 
outer horizon the field is shown as a function of v (bottom curve). For the clarity of the figure 
we display only the maximas of the oscillations (the field oscillates with a period of T = ir/m to 
within 0.2%). At late-times the field's amplitude decays slower than any power-law. The initial 
data are those of Fig. |[ 
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